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Abstract 



We show that QCD observables defined as infrared- and collinear-safe, according 
to the usual Sterman- Weinberg criteria, can nevertheless be infinite at accessible 
points inside phase space, to any finite order of perturbation theory. Infrared 
finiteness is restored after resummation of divergent terms to all orders. The 
resulting characteristic structure, which we call a Sudakov shoulder, represents 
an interesting new class of QCD predictions. 
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1 Introduction 



The perturbative QCD approach P] to the calculation of hadronic cross-sections at large 
momentum transfer Q is based on the Sterman- Weinberg criteria || of infrared and collinear 
safety, combined with the factorization theorem for mass singularities A hadronic 
observable is said to be infrared and collinear safe if it is insensitive to the emission of 
soft momenta and to the splitting of a final-state momentum into collinear momenta. The 
properties of infrared and collinear safety guarantee that the long-distance component of 
the scattering process, which is controlled by non-perturbative phenomena, is suppressed by 
some inverse power of Q as the hard scale Q increases. Thus, the observable is dominated 
by short distances and the short-distance component, which depends only logarithmically 
on Q, is computable as a power-series expansion in the strong coupling as(Q)- 

It is commonly assumed that the consequences of the Sterman- Weinberg criteria are 
even wider: the power-series in as is assumed to have coefficients that are finite order by 
order in perturbation theory. Interpreted in this way, the criteria become a statement on 
the infrared and collinear finiteness of the fixed-order perturbative expansion. 

In this paper we would like to point out that, in general, this strong interpretation of 
the criteria is not valid. In fact, the perturbative expansions of infrared- and collinear-safe 
quantities can be divergent order by order in perturbation theory and these divergences are 
still produced by the radiation of soft and/or collinear partons. 

The presence of logarithmic singularities that spoil the convergence of the perturbative 
expansion near the exclusive boundaries of the phase space of infrared- and collinear-safe 
observables is well known These singularities reflect the difficulty in extending the 



Sterman- Weinberg criteria towards extreme kinematical regions and, in many cases, they 
can be handled by all-order resummation methods. 

The main issue to be discussed in the present paper concerns a new general class of 
singularities that appear inside the physical region of the phase space. They arise whenever 
the observable in question has a non-smooth behaviour in some order of perturbation theory 
at an accessible point, which we shall call a critical point, inside phase space. This can 
happen if the phase-space boundary for a certain number of partons lies inside that for a 
larger number, or if the observable itself is defined in a non-smooth way. In either case, if 
the distribution of the observable is actually discontinuous at that point in some order, it 
will become infinite there to all (finite) higher orders. If it is continuous but not smooth, 
i.e. if some derivatives are discontinuous, then the higher-order predictions will become 
unstable, with infinite derivatives, at that point. 

There are several possible reactions to this state of affairs. One can try to avoid non- 
smooth observables altogether. As we shall see, this actually rules out some commonly-used 
quantities. One can use such observables but avoid the phase-space regions near critical 
points. In this case, one still needs some understanding of the phenomenon in order to 
assess the extent of the dangerous regions. Alternatively, one can identify the problematic 
terms in each order of perturbation theory and resum them to all orders. As mentioned 
above, this approach has been used successfully for the treatment of singularities near the 
exclusive phase-space boundary in many observables. We advocate a similar approach to 
the new class of singularities discussed here, although much work remains to be done in 
order to reach a comparable level of understanding. As a first step in this direction, we 
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carry out the resummation of leading double-logarithmic terms, and show that it leads to 
smooth behaviour at the critical point. 

In the following Section we give a general discussion of the effects of soft and/or collinear 
parton emission in fixed-order QCD perturbation theory. We first review the cancellation 
mechanism which normally ensures the regularity of infrared- and collinear-safe observables, 
and the way in which this breaks down at the exclusive boundary of phase space. We then 
show how a similar breakdown will occur inside phase space when the observable is defined 
in a non-smooth way, and describe the general form of the resulting singularities. 

In Sect. |3]we discuss in some detail two specific cases in which such singularities actually 
appear. The first is the distribution of the C-parameter [T2],I3], a well-known event shape 
variable for e + e~ annihilation final states. The physical region is < C < 1, but at the 
point C = 3/4 there is a divergence in the perturbative prediction of the C-distribution 
at second and higher orders. The second example is taken from jet physics in hadron- 
hadron collisions, and concerns the jet shape or profile function, which describes the angular 
distribution of energy with respect to the jet axis [14]. Here one finds a critical point at 



the (arbitrary) boundary of the cone used to define the jets, which leads to divergences in 
jet shapes and exclusive multijet cross sections at next-to-leading order and beyond. 

In Sect. |] we discuss the resummation of these divergences to all orders. We expect on 
general grounds that infrared finiteness will be restored by resummation, and we show that 
this is the case for the leading double-logarithmic terms. After resummation, one obtains 
a characteristic structure which is not only finite but smooth (infinitely differentiable) at 
the critical point. By analogy with the characteristic Sudakov peak generated by a similar 
mechanism near the exclusive phase space boundary, we call this structure a Sudakov 
shoulder. We also discuss the resummation of a more general class of non-smooth behaviours 
at the critical point, which give rise to different structures related to the Sudakov shoulder. 

Finally in Sect. [5] we summarize our main results and conclusions. 



2 Singularities in fixed order 

The singularities we are interested in are due to the radiation of soft and/or collinear 
partons. In perturbative calculations, the soft and collinear regions of the phase space 
lead to infrared divergences in the QCD matrix elements. In the case of infrared and 
collinear safe observables the divergences cancel upon adding real and virtual contributions 
but, precisely speaking, the cancellation does not necessarily take place order by order 
in perturbation theory. In fact, in particular kinematic configurations real and virtual 
contributions can be highly unbalanced, spoiling the cancellation mechanism. 

Soft-gluon cancellation 

In order to discuss this point most generally, we denote by C an infrared- and collinear- 
safe observable defined in a hard-scattering process at the scale Q. We can always assume 
that C is dimensionless and varies in the kinematic range < C < 1. We denote by o^(C) 
a related cross section or distribution^ evaluated at ra-th order in perturbation theory and 

^The perturbative distribution a(C) can be either a measurable quantity or a short-distance coefficient 
function that appears in the factorization theorem. 
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consider the effect of soft-gluon radiation at the next perturbative order. If 1 — z denotes 
the fraction of the centre-of-mass energy involved in the radiative process, virtual (v) and 
real (r) soft gluons affect the cross section with the following emission probabilities 

dw v (z) _ a s (Q) 5{1 _ z) ( X -*J?L ln ^_^ {21) 



dz ix Jo 1 — z' 1 — z 

= ^MJ_i n -J_e(i-,- e ). (2.2) 

dz 7T 1 — z 1 — z 



The expressions ( 0| , |2.2j) are valid to double-logarithmic accuracy and arise from the com- 
bination of the customary bremsstrahlung spectrum du/uo with the angular distribution 
d8 2 /8 2 for collinear radiation. Here we have introduced an unphysical cutoff e on the min- 
imum energy fraction of both gluons because the probabilities in Eqs. ( |2.1|) and (|2.2j ) are 
separately divergent if the cutoff is absent. However, in the case of infrared- and collinear- 
safe observables, real and virtual emissions contribute with equal weight K to cross sections, 
although the actual value of the coefficient K depends on the particular observable under 
studyQ Thus, adding the real and virtual terms, the physical limit e — > can be safely 
taken, leading to a finite differential probability: 



dz e-o 



dw v (z) dw r (z) 



(2.3) 



dz dz 
In-!-) , (2.4) 



X-z 

where a = KasiQ)/^ an d, as usual, the plus-prescription [g(z)] + stands for 

1 dzf(z) [g(z)] + = f dz [f(z) - /(l)] g(z) . (2.5) 
Jo 



The plus-prescription in Eq. fl2.4j) is the result of the cancellation of the dependence 
on the unphysical cutoff e and defines a well-behaved distribution, i.e. one leading to non- 
singular functions, when acting on any smooth function f(z) that is finite at z = 1. How- 
ever, this is the sole cancellation mechanism of infrared singularities that is guaranteed by 
infrared and collinear safety in fixed-order perturbation theory. The Sterman- Weinberg 
criteria do not imply that the coefficients of the perturbative expansion are non-singular 
functions. In fact, normally these coefficients are (or are obtained from) singular generalized 
functions or distributions, like the soft-gluon probability in Eq. fl2.4|) . Such distributions 
lead to finite quantities only when they are integrated with sufficiently smooth test func- 
tions. 

In particular, the (n + l)-st order cross section a^ n+l \C) is obtained by combining the 
soft-gluon probability with the n-th order term a^ n \C) and the latter may well not be a 
smooth function of C for all values of C . In this case, the real and virtual contributions are 
effectively unbalanced and singularities show up in the fixed-order perturbative expansion. 

In order to discuss these singularities, we have to specify how the observable C depends 
on the soft-gluon energy. We consider for definiteness the case in which C increases when 



■'"In the case of measurable quantities, the coefficient K is typically positive definite. If instead the 
distribution <r(C) is a short-distance coefficient function, the sign of K depends strongly on the factorization 
scheme used to define the parton distribution functions. 
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soft partons are radiated in the final state. Thus, the emission of a soft gluon with energy 
fraction 1 — z produces a change C —* C + 5C, where 5C — 0(1 — z), and, without loss of 
generality, we can write the next perturbative order as follows: 



a 



^{C) = f dya^iC -y) (*^) +... , (2.6) 

" V / 2 = 1 — J/ 



where the dots stand for less singular terms, if any. Because of the character of the soft- 
gluon distribution dw(z)/dz, if is not smooth at some point in the kinematic range 
< C < 1, the integral on the right-hand side of Eq. (|2.6|) can be divergent. 



Singularities at the exclusive boundary 

The most common and extensively studied soft-gluon singularities are those arising at 
an exclusive boundary of the phase space. In this kinematic regime, the radiative tail of real 
emission is strongly suppressed, producing the loss of balance with the virtual contribution. 
Since the observable C that we are considering increases when soft partons are radiated, 
the exclusive boundary is the region C — > 0. Here the non-smooth behaviour of a^ n+1 ' ) (C) is 
simply due to the associated phase-space constraint 6(C). Inserting Eq. (|2.4|) into Eq. 
we have: 



a 



("+!)(C) = a 



--a^(C) \n 2 C+ / C ^ln- (a^(C - y) - a^(C)) 
2 Jo v v ' 



+ ... . (2.7) 



Thus, even if a^ n \C) is finite^ as C — > 0, the first term in the square bracket produces a 
divergent contribution at the kinematical boundary: 

o^\C) c ~ o -|aW(C)ln 2 C + ... . (2.8) 

These double-logarithmic terms are usually called soft-gluon singularities of Sudakov 
type. They appear in the perturbative expansions of many e + e~ shape variables in the 
two-jet limit ||||. The observable C, for instance, can be the C-parameter [|T2"| , |13 



or 



C = 1 — T, where T is the thrust |T5[ . Similar soft-gluon effects occur in hadron collisions 



for Q^-distributions in the Drell-Yan process [C ~ Q±/Q) and for the production of 
systems of high mass M near threshold (1 — C ~ M/y/S, where \/~S is the centre-of- 
mass energy). Outstanding examples of these systems are lepton pairs with large invariant 
mass produced via the Drell-Yan mechanism 0, the hadronic final state in deep- inelastic 
lepton-hadron scattering 0, heavy quark- ant iquark pairs [P-fTTH and pairs of jets at large 



transverse momentum 11 



Singularities inside the physical region 

Another possible source of soft-gluon singularities, which has received less attention in 
the literature, arises when a^ n '[C) has non-smooth behaviour at a certain value C = Co 
inside the physical region of phase space. This can happen if the phase-space boundary for 
a certain number of partons lies inside that for a larger number, or if the observable C is 
defined in a non-smooth way. We discuss examples of both types below. 

§If <j( n )(C) is divergent at C — 0, the double-logarithmic singularity in Eq. ( pTg| ) is enhanced and can 
also be non-integrable. 
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The diagnosis of singularities inside the physical region is physically less obvious than 
that for singularities at the phase-space boundary. At the exclusive boundary we are dealing 
with a multi-scale process: besides the hard scale Q, there is also another natural scale, 
the inelasticity scale Q' ~ CQ, that plays a relevant role. The Sudakov singularities in 
Eq. (|2.8|) follow from the large mismatch between these two scales, Q' <^ Q when C — > 0. 
For singularities inside the physical region, on the other hand, the lack of balance between 
real and virtual contributions is just produced by the sharpness of the distribution around 
the critical point C , and the actual identification of a critical point requires careful analysis 
of the kinematics and dynamics. 

To show how these singularities arise, we denote the cross section for C < Co by <7_(C) 
and that for C > C by a + (C), and we can suppose that at the n-th order both and 
<7+^ are infinitely differentiable functions. Then we find from Eq. ( |2.7| ) that a^ +l \C) is 
infinitely differentiable and regular at C = Cq while a+ (C) can be singular there. The 
difference comes from the second term in the square bracket of Eq. ( p.7|) and is given by 

a ( : +1) (C) - a^ +1 \Co) = -\ ln 2 (C - C ) { [<rf(C7 ) - at\C ) 

+ (C-C Q )[a^'(C )-a^'{C )]+...}+... , (2.9) 

where a^'(C) denotes the first derivative of the distribution with respect to C and the dots 
indicate terms that are less singular as C — > C . Thus if a+ (Co) ^ a_ (Cq), i.e. a^(C) 
has a step at C = Co, then a+ {O) has a double-logarithmic divergence at that point. 
The same mechanism leading to Eq. (|2.9| ) will enhance the double-logarithmic singularity 
by further integer powers of ln(C — Co) in yet higher orders^]. 

In general, a step at some point C = Co inside the physical region at the n-th order of 
perturbation theory always generates more and more divergent contributions at C = Co in 
higher and higher orders. To any fixed order m > n, the cross section a^ 1 (Co) is infinite, 
because the emission of arbitrarily soft gluons is not cancelled by virtual corrections at that 
point. One may regard this as a failure of the Sterman- Weinberg criteria for infrared and 
collinear finiteness in fixed-order perturbation theory: we have considered an observable 
which satisfies the criteria, but nevertheless the perturbative prediction at a point inside 
the physical phase-space is divergent. 

One might think that the construction of an infrared and collinear safe quantity, which 
at the n-th perturbative order has a step-like behaviour inside the physical region, is quite 
abstract. In fact, this can happen if, due to the finite number of partons involved in the 
calculation at that order, they can occupy the phase-space region only up to a value C = Co 
that is smaller than the maximum value permitted by kinematics. To produce a step it is 
sufficient that the computed distribution does not vanish at C = Co- As we shall see in 
Sect. ||, the distribution of the C-parameter in e + e~ annihilation has these features. 

A similar behaviour can occur if the infrared and collinear safe observable C is used 
to classify the topology of hadronic events. In this case sharp structures can easily be 
produced. For instance, C can be a jet resolution parameter that is defined in a non- 
smooth way (e.g. it can have a non-trivial and non-smooth dependence on energies and 
angles of partons). Then the topology of the (partonic) event, e.g. the number of jets, 

^Note that all these singularities are integrable in a neighbourhood of the critical point. 
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can suddenly change at a certain value Cq. At that point the rate a^ n \C) for producing 
a fixed number of jets can have a step. Moreover, a corresponding step, with the opposite 
discontinuity, will appear in the production rate relating to a different number of jets. At 
the next perturbative order, the step in a^ n \C) leads to a singularity for C > Cq and, 
by the same argument as in Eq. (|2.9|) , the step in the other jet rate leads to a similar 
singularity for C < Cq. Since the sum of the two jet rates remains finite, both of them 
will show a logarithmic divergence on both sides of the step, that is, above and below the 
critical point C = Co- This may appear in contradiction with the statement above Eq. Q2.9|) . 
However, since we are considering a variable C that is defined in a non-smooth way, it is not 
monotonic with respect to soft-gluon energies. Therefore, the simple formula in Eq. (|2.6| ) 
cannot be applied, although the general relation between steps and singularities is still 
valid. A phenomenologically-relevant example of this kind of singularity is also considered 
in Sect. 

The result in Eq. Q2.9Q also allows us to discuss the effect of soft-gluon radiation on non- 
smooth quantities when they do not have a step-like behaviour. For example, the n-th-order 
distribution cr^(C) may be continuous but not its first derivative (a+ (Co) 7^ <r_ '(Cq)), 
so that a^{C) has a sharp edge at C = Co. An example would be C = 1 — T where T is 
the e + e~ thrust variable, which has kinematic range \ < T < 1 but is zero in first order 
(n = 1) for T < I and vanishes linearly as T — > To = | + . Then a^ n+1 \C) is also continuous 
but has a cusp at C = Co, where a+ (C) diverges double-logarithmically[|. In the case 
of C = 1 — T, from Eq. Q2.SQ we obtain that a^'(C) diverges to negative infinity^ at 
Cq = 1 — To — 3- 111 higher orders <r(C) remains continuous but the degree of divergence 
of its derivative is enhanced by an increasing number of double logs, so the edge gets more 
and more like a step order by order in perturbation theory. 

An edge at a certain order thus does not produce divergences in the perturbative expan- 
sion. Nonetheless, fixed-order predictions are likely to be unreliable near the critical point 
because the shape of the distribution, independently of its magnitude, is highly unstable 
with respect to radiative corrections. For instance, if a^ n \C) is decreasing at the edge and 
the perturbative coefficient a in Eq. ( |2.9| ) is positive, the cusp at the next order leads to an 
increasing distribution for C > C . 

3 Examples 

Here we illustrate the above general discussion with two specific examples of distributions 
of well-known observables which do develop divergences at points inside the physical phase 
space, due to imperfect cancellation between real and virtual soft-gluon singularities. 

3.1 The C-parameter in e + e~ annihilation 

The C-parameter is an infrared- and collinear-safe observable that was introduced in 
Refs. JT^] to describe the shape of hadronic events produced by e + e~ annihilation. De- 

"in general, if is continuous up to its fc-th derivative, the double-logarithmic singularity affects the 
fc-th derivative of a^ n+1 \ 

**This cusp is visible in Fig. 1 of Ref. @. 
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noting by pf the components of the three-momentum of any final-state particle i, one first 
considers the linearized momentum tensor 

r\a/3 _ J2iP?Pi '' /\Pi\ /o i \ 

and computes its eigenvalues Aj. The C-parameter is then defined as follows 

C = 3(A 1 A 2 + A 2 A 3 + A 3 A 1 ) . (3.2) 
By definition the eigenvalues satisfy the following constraints 

< A, < 1 , £ Ai = 1 , (3.3) 

i 

so that the C-parameter varies in the kinematic range < C < 1. In particular, C = for 
a perfectly two-jet-like final state (e.g. A 2 = A 3 = 0, Ai = 1) and C — 1 for an isotropic 
and acoplanar distribution of final-state momenta (Ai = A 2 = A3 = |). 

In fact, the maximal value C = 1 can only be achieved when there are four or more 
final-state particles. Planar events have one vanishing eigenvalue and occupy the kinematic 
region C < |. Hence, we are in a situation where the phase-space boundary for three 
particles is below the kinematic limit of the shape variable. For a general three-particle 

state 

C = C 3 (xi,x 2 ) = 6(1 -xi)(l -x 2 )(l - x 3 )/{x x x 2 x 3 ) (3.4) 

where x« = 2pi ■ Q/Q 2 are the centre-of-mass energy fractions (xi + x 2 + x 3 = 2). The 
maximum value C — | for three particles corresponds to the symmetric configuration 
Xi = x 2 = x 3 = §. 

In perturbation theory, the distribution of the C-parameter for C ^ has the general 
form 

-% = a s A(C) + a 2 s B(C) + 0{a%) , (3.5) 
where a$ = as(Q)/27i and we normalize to the Born cross section a , as was done in 



Ref. |T7|]. The two functions A(C) and B{C) can be identified respectively with the dis- 
tributions and a^ n+1 ^ of Sect. ||]. In particular, we are exactly in the situation where 
Eq. ( |2.6| ) applies because, as can be checked easily, the C-parameter increases when addi- 
tional soft particles are produced. 

The first-order distribution A(C) is given by the qqg final state: 

= / dx\dx 2 Q(xi + x 2 — 1) Mixx, x 2 ) S[C — C 3 (x\, x 2 )\ (3.6) 
Jo 

where 

M {xuX2 )=C F (1 _ X l^_ x2) . (3.7) 

Note that the matrix element ( |3.7| ) has a finite value M = M(x% = x 2 = |) = 8Cp in 
the symmetric configuration where C — |. Furthermore one finds that the available phase- 
space remains finite as C — > | . Therefore, the 0(as) distribution ( |3.6|) has a stepQ at 

^ Note that, consistently with the general discussion in Sect. ||, the first-order distribution A(C) in 



Eq. (^6) diverges at the exclusive boundary C = of the phase space. This divergence and the corre- 



sponding ones in higher orders can be resummed |f3| using the techniques of Ref. 
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Figure 1: Predictions of the C-parameter distribution for as 
Solid: 0(a s ). Dot-dashed: resummed. 



0.12. Dashed: 0{a s ). 



the three-parton upper limit C 



3. 
4' 



^(C) ~ 3 A(f)e(f-C) , (3.8) 

O ~ -J/4 

A(f) = 243^^ = 243^^. (3.9) 

3" 



The derivative of A{C) is also finite at C — | : in fact 

^(!) = -H!) • (3-io) 

The resulting first-order prediction in the vicinity of C = | is shown by the dashed curve 
in Fig. 

According to Sect. 0, the second-order contribution -B(C) should be smooth as C — > § 
but should have a double-logarithmic singularity as C — > | + . This is indeed the case, as 
shown by the solid curve in Fig. [I]. Denoting -B(C) above/below C = f by B±(C) we find 
using the Monte Carlo matrix element evaluation program EVENT [[HJ 

= 230 ±10 . (3.11) 

Performing the analytic calculation for C — > | + , we find 

B+(C)~A(|) [(2C F + C A )(l-5)ln 2 5+(3C F + ^C A -|Ar / )ln(5] +/i(C) (3.12) 

where 

*=§(C-f) (3.13) 
and the remainder is expected to take the form 

/i(C) = /i(!) + 0(<J1ii<$). (3.14) 
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L = -log(8C/3 - 2) 
Figure 2: Second-order prediction of the C-parameter distribution for C 



Points: EVENT Monte Carlo. Curve: Eqs. (gTT^-gTTS]) 



Using Eq. (|3.10 ), we see that the first term in the square bracket on the right-hand side 
of Eq. fl3.12j ) is in complete agreement with Eq. ( [2.9|) with a = (2CV + Ca)ois{Q)/t{. The 
second term in the square bracket is due to collinear but non-soft parton splitting: single- 
logarithmic contributions of this type were neglected in the simplified discussion of Sect. 0, 
which was valid to double-logarithmic (DL) accuracy only. 

Comparing our analytic calculation with numerical data from EVENT |[L7| , we find 
good agreement, as shown in Fig. §. The value of the constant /i(|) in Eq. ( |3.14D is found 
to be 

h{\) = 146 ±3 . (3.15) 



3.2 Jet shapes in hadron collisions 

As our second example we study the angular distribution of transverse energy (Et) in jets 
produced at high Et in hadron-hadron collisions ||14|| . Here it is customary to define a jet 



by a cone algorithm ||18|| , that is, by maximizing Et with respect to the direction of a cone 
of opening angle R in pseudorapidity, azimuth (77, </>) space. The cone axis then defines 
the jet axis (rjj, < 



!>j). Alternatively a /c_|_-clustering algorithm can be used |19[]. In this case 
again a parameter R can be introduced, effectively representing the jet cone size 



Several different versions of the cone algorithm are actually used both in theoretical 
calculations and in experimental analyses. They differ with regard to minor details and 
sometimes major points. The latter are related to infrared and collinear safety. In some 
versions [^TJ, the direction of the unassigned cluster (parton) with the highest Et is used 
to start the clustering procedure. Since the highest-Er parton is not stable with respect to 
collinear splitting, these algorithms are collinear unsafe. Their unsafeness certainly affects 
any perturbative calculation that involves more than 5 partons in the final state. Collinear 
safety is fulfilled by iterative-cone algorithms that use all clusters, possibly excluding those 
with E T below a fixed threshold value [P^] , as seed directions of the cone. Nonetheless, these 
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algorithms are still infrared unsafe in the absence of an E^-threshold [23|, and therefore 
strongly affected by the threshold when one is applied [p4" |. In iterative-cone algorithms, 
infrared and collinear safety can be achieved [^5|] by using the ^-weighted midpoints of all 
pairs of jets (partons) as additional seed directions. In the following, by 'cone algorithm' 
we refer to these improved versions. 

We consider here the differential jet shape or profile function ip(r), defined such that 
ip(r) dr is the fraction of transverse energy lying within a ring of radius r, width dr, centred 
on the jet axis. We have 

,/ s Ei E Ti S(r-Rjj) 

W = tp Q / D > ( 3 - 16 ) 



J2k E Tk Q{R — R 



kj) 



where 



R u = vim ~ Vj) 2 + (4>% - 



(3.17) 



and the sums over i, k can be over either all particles or only those particles assigned to 
the jet. The function ip(r) has been studied in detail in Ref. [24], and many of the points 
we make are also discussed there. 

For r > 0, if)(r) starts at order as and at this order its r-dependence is due to 3-parton 
final-states. Let us simplify the discussion by considering a cone size R < tt/3 and three 
partons of momenta Pi,P2,P3 with transverse-energy ordering, e.g. En > E?2 > £t3- In 
this case there are two possible jet configurations: 

{1},{2},{3} 



a) jets 



b) jets = {1}, {23} 

yalg 



if Rf 3 9 > R 
if R a 2 l 3 9 <R 



where i? 23 depends on the jet-algorithm: 



ricone 
K 23 



E 



T2 



Et2 + E< 



-R 



2.S 



R 



T3 



23 



R 



2?> 



(3.18) 



The jet shape ip{r) can be written in terms of two contributions ipi n {r) and ipouti^), arising 
from particles assigned and not assigned to the jet, respectively. To order as, the two 
contributions are proportional to the following expressions 



oc e(R-R$)6(J- {23}) 

Et3 



E T2 8 ( r 



Et2 + E> 



R 



2.3 



T3 



+ E T3 5 f r 



E- 



T2 



Et2 + E'_ 



-R 



2.3 



T3 



tfi(r) oc Q{Rfi - R) [5{J - {2}) E T3 + 6(J - {3}) E T2 ) S(r - R 23 ) 
where we use the notation 



(3.19) 
(3.20) 



5{J - {2}) = 5{E TJ - E T2 ) S(rjj - rfy) 5{<j> 



(3.21) 



and 
have 



and likewise for 8{J — {3}) and 5(J — {23}). Equations ( |3.19| ) and ( |3.20| ) show that ip^ (r) 



are respectively proportional to Q(R — r) and 0(r — R), and when r — > R we 



4l\R)>0, ^i(r) 







(cone alg. 



out ( 



VffiOR) > (k ± alg. 



(3.22) 
(3.23) 
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In conclusion, i[)^(r) in the cone algorithm has a step at r = R, while in the k± algorithm 
it has an edge (step) if only particles assigned to the jet (all particles) are considered. 

Let us now investigate possible singularities at the next order. Since ip^ 2 '(r) oc -Eri 
S(r — Ru), we should consider how the angular distance between the jet direction and 
that of parton i is affected by collinear splitting and soft gluon radiation. In the case 
of collinear splitting the angular distance does not change. Thus we have to examine 
soft-gluon radiation with energy fraction 1 — z and single-logarithmic emission probability 



dw(z) 
dz 



In the cone algorithm, Rij is always smaller than R if i belongs to the jet and always 
larger than R if it does not. If i belongs to the jet and a soft gluon is emitted inside the 
cone, cancellation of real and virtual divergences will occur as usual, leading to a smooth 

(2) 

contribution to ij)\ n . However, if the gluon is emitted outside the cone there is a mismatch 
of real and virtual contributions, which will generate a divergence at the boundary of the 
cone. The corresponding change in R^j is — c(l— z), where c > depends on the kinematics. 
Similarly, if parton i is outside the jet and a gluon is emitted into the cone, the shift in Rij 
is +c(l — z). Thus, according to our general discussion in Sect. ||, we find 

4*\r) ~ Q(R-r) {-a 4$(R) In + const.} 

^S(r) ~ 6(r - R) const. (3.24) 

The coefficient a is likely to be positive in most jet configurations so that i/j( 2 \r) diverges 
to — oo when r — > R from below and stays finite when r —> R from above. 

In the k± algorithm at higher orders, because of the iterative recombination procedure, 
Rij can be either larger or smaller than R, independently of whether i is or is not a particle 
in the jet. Thus both and Vowk r ) are non- vanishing in a neighbourhood of r = R. 

Moreover, Rij does not vary monotonically with z when a soft gluon is added to the jet. 
Thus we find that ipl n (r) has a cusp at r — R and ipout( r ) nas the following double-sided 
singular behaviour^ 

^(r) ~ a \q(R - r) ^ t (R) In — - 6(r - R) 4Ji(R) In + const.) . (3.25) 

I R — r r — R J 

The double-sided singularity follows from our general discussion in Sect. |^ on non-monotonic 
observables with stepwise behaviour at lower order. 



Note that the singularities in Eqs. ( |3.24| ) and ( |3.25| ) are both single-logarithmic. This 



is because they are due to soft but non-collinear radiation. 

Related problems connected with the calculation of jet cross sections in higher orders 
have been discussed in Ref. |[23|| . The definitions of the fixed-cone and iterative-cone al- 
gorithms given there are actually unsafe beyond a certain order. More precisely, since 
the highest unassigned .Sr-cluster (parton) is used to start the clustering procedure, these 
algorithms are collin ear unsafe. For a cone size R < 7r/3, their unsafeness affects perturba- 
tive calculations that involves more than 5 final-state partons. Thus one may still perform 
calculations up to order a| using these definitions, provided that divergences at critical 
points are carefully handled. 



Assuming a > 0, this agrees with the results of Ref. pi 
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Consider the 3-jet cross section in the iterative cone algorithm as defined in Ref. [[23] at 
a fixed value of R23, the angular distance between the two less energetic jets. At LO this 
cross section has a step at the critical point R23 = R, the cone radius. Actually it vanishes 
for R23 < R and is finite and not vanishing for R23 > R. At NLO the emission of a soft 
gluon at a large angle produces a single-logarithmic divergence (the soft gluon not being 
collinear) at -R23 = R, of the type discussed above. 

This divergence is integrable and it does not produce any singularity in the cross section 
integrated over R23. On the other hand, if the phase-space slicing method |26| is used to 



perform the calculation, the integrable divergence can give rise to a lns min dependence, 
where s m i n is the 'slicing parameter' used for the numerical evaluation of the cross section. 
Such a dependence can be numerically important, even though it should disappear for 
sufficiently small values of s m i n . 

To illustrate this point we note that the cutoff e introduced in Sect. |2| to define the 
plus-distribution is somewhat analogous to a slicing parameter. In the case of a stepwise 
behaviour of fr( n > at C = Co, the integral up to C\ > Co, 



lim [ * dC fr( n+1 \C;e) , (3.26) 

e—0 J C 



is in fact finite, but if we keep e fixed we get for C\ — Co < e 



[ Cl dC a {n+1) (C;e) = a -^ln 2 e / * dCa { +\C) + const. + 0{e In e) 

J Co ^ J Cn 



I C 



(3.27) 



which shows an unphysical dependence on lne. The unphysical logarithmic dependence 
does not cancel as long as e is not sufficiently small. 



4 Resummation and the Sudakov shoulder 

Although we have shown that divergent fixed-order predictions will generally arise at non- 
smooth points inside the physical region, one might feel that this is a purely academic 
difficulty, since the divergence at such a point C = C is integrable, as long as a^ n \C) 
itself is integrable in a neighbourhood of that point. The cancellation of real and virtual 
divergences will therefore take place in the presence of smearing in C, as would be expected 
from non-perturbative (hadronization) effects. However, the resulting prediction will be 
highly unstable with respect to the order at which we stop the perturbative calculation 
and the amount of smearing introduced. 

In particular, since the hadronization smearing should cancel divergent terms pro- 
portional to some power of as{Q), this would imply that non-perturbative effects scale 
logarithmically with Q, thus spoiling not only the finiteness but also the safety of the 
Sterman- Weinberg criteria. We should then be forced to conclude that some infrared- and 
collinear-safe observables are affected by non-perturbative contributions that are not power 
suppressed. 

Alternatively, one might think that the divergence in the fixed-order prediction can be 
cancelled by the resummation of still higher-order perturbative corrections. One might 
argue that the resummation of real and virtual terms produces a Sudakov suppression of 
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the divergence, similarly to the Sudakov suppression of many observables at the exclusive 
boundary of the phase space This scenario would lead to a resummed prediction that, 
at the critical point C = Co, is finite but still has a stepwise behaviour. On physical grounds 
we do not anticipate this behaviour, and we should expect non-perturbative effects to fill 
up the gap between the two sides of the step. Since the gap is produced in perturbation 
theory, in this case also the necessary non-perturbative contributions could not be power 
suppressed. 

We believe that instead the problem of divergences inside the physical region has a 
satisfactory solution entirely within the context of perturbation theory. Namely, the re- 
summation of the soft-gluon contributions to all orders, rather than producing a Sudakov 
suppression of the divergences, leads to a structure that is continuous and smooth, indeed 
infinitely differentiable, at C = Cq. We call this structure a Sudakov shoulder. The ex- 
istence of a Sudakov shoulder implies the restoration of validity of the Sterman- Weinberg 
criteria at infinite order in perturbation theory. 

To demonstrate this point we first recall how soft-gluon resummation produces the 
Sudakov form factor at the exclusive boundary of the phase space. We thus consider the 
case in which cr(C) is a measurable quantity (with a positive soft-gluon coefficient a) and 
proceed to the iterative application of Eq. ( |2.7|) , retaining for simplicity only the double- 
logarithmic (DL) contribution at each order. Taking into account the symmetrization with 
respect to exchange of soft gluons, we obtain an exponential series that, after summation 
in the neighborhood of the point C = 0, gives 



a 



(00) 



(C) c = ex P {-|ln 2 C 



a 



(no) 



(C) 



(4.1; 



where hq represents the lowest perturbative order for our observable. The exponential term 
on the right-hand side of Eq. (|4.1|) is the customary Sudakov form factor, which suppresses 
the observable at the exclusive boundary. 

The analogous iterative procedure in the region C ~ Co also leads to an exponential 
series which can be summed to give 



a 



(oo) 



(C) 



a 



(oo) 



(Co) 



expj-^ ln 2 (C-C )} 



a 



(™o) 



(C 



a 



Oo) 



(Co) 



(4.2) 



where Uq represents the order at which the step at C = Co first appears. Now we may 
assume that cr*L°\Co) has no soft-gluon singularities, since by hypothesis C is always in- 
creased by soft gluon emission. It then follows from Eq. (|4.2|) that after resummation the 
cross section is finite, continuous and infinitely differentiable at C = Co. Rather than sup- 
pressing the divergence at C > Co, Sudakov resummation leads to the suppression of the 
step. The general form of the resummed cross section will be a smooth extrapolation from 

r (oo) 
r a (oo) 

of cri 



rr _ (Cq) into the region C > Co, joining smoothly with the different asymptotic behaviour 



(C) 
shoulder. 



well above C = Co. This is the characteristic structure we call a Sudakov 



The dot-dashed curve in Fig. [I] illustrates the Sudakov shoulder in the case of the C- 
parameter distribution. Here we simplify the treatment of the region C > § not only by 
using the DL approximation for the form factor but also by approximating the lowest-order 
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contribution A(C) in Eq. (|3.6|) by a step function, A(C) ~ A(| )©(§ — C). We then obtain 



1^ ~ a 5 A(|) {9(1 - CO + e(C - 1) (l - exp [-2A«a 5 ln 2 (C - §)] ) } (DL) 



(4.3) 

= C f + \C a . (4.4) 



adC 
where 



Expanding the right-hand side of Eq. ( [4. 3D perturbatively, one finds a series of logarithmi- 
cally divergent (although integrable) terms for C — > f + . The summation of these terms to 
all orders leads to a finite result and the lowest-order step-like behaviour is smoothed. The 
resummed expression (|4.3|) shows a shoulder that extends beyond the value C = |. The 
shoulder becomes less steep with increasing or as- Since = Cp + Ca/2 ~ 2Cp, 
in an abelian theory (i.e., setting Ca = 0) the DL shoulder would be twice as steep as in 
QCD. The additional broadening in the QCD case is due to the fragmentation of the gluon 
jet. Note that the shoulder is nevertheless still quite sharp on the scale shown in Fig. 1. 

Soft-gluon resummation will smooth out perturbative predictions of infrared- and collinear- 
safe observables at other types of critical points inside phase space. After the case of a 
step, the next simplest example is that of a critical point C = Co with an edge at a given 
perturbative order no- The resummed behaviour at this point is easily obtained by noting 
that the first derivative &'(C) of the corresponding distribution is a safe quantity that has 
precisely a step at order no- Thus we can apply to a'{C) a DL resummation analogous to 
that carried out in Eq. (|4.2|), to obtain the following distribution 



a^°\C) = a { f°\C ) + (C - C ) ( xi oo), (Co) + \a^ o) '(C ) - ai no)/ (C )l \° °° dxe~^ x . 



(4-5) 

The first two terms on the right-hand side extrapolate <r_ (C) linearly into the region 
above the critical point C = Co and the last term bends this linear extrapolation smoothly, 
in proportion to the sharpness of the lower-order edge. 

The resummation of double-sided singularities of the type encountered in jet profiles is 
technically more complicated. It leads to a smooth "jump" structure obtained by smearing 
(convoluting) two Sudakov shoulders: a shoulder above and an inverted shoulder below the 
critical point, or vice versa. As a result, the two sides of the lower-order step match with 
some intermediate value at the critical point in the all-order distribution and the step is 
smoothed. An example of this resummation carried out numerically for the jet shape ip(r) 
is presented in Ref. p4j (cf. Fig. 11). In general, the derivation of resummed expressions 
in analytic form is quite difficult and not always feasible. This is because the type of 
smearing (convolution) to be applied to the two shoulders is strongly dependent on the 
detailed kinematics. In the presence of double-sided singularities, the kinematics is usually 
very involved since it is responsible for the non-monotonic behaviour of the observable with 
respect to variations of particle angles and energies. 



5 Discussion and conclusions 

In this paper we have discussed in rather general terms the divergences which can arise in- 
side the physical region in infrared- and collinear-safe quantities computed to any fixed order 
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in perturbation theory. We have also given specific examples of widely-used quantities in 
e + e~ and hadron-hadron physics for which these problems actually arise. The divergences 
correspond to integrable singularities and therefore they could in principle be removed by 
non-perturbative smearing effects. However, this would require non-perturbative contribu- 
tions that are not power-suppressed at high energies. 

We have argued that instead the remedy for these problems lies entirely within perturba- 
tion theory: the resummation of soft-gluon contributions to all orders should be sufficient to 
render any safe quantity finite and smooth throughout the physical phase space. Despite 
the violation of finiteness of fixed-order calculations, perturbative resummation suggests 
that non-perturbative effects are still power-suppressed in infrared- and collinear-safe ob- 
servables. 

A resummation of double-logarithmic terms was performed to illustrates how this works 
in the case of the e + e~ C-parameter distribution (Fig. 1). More work remains to be 
done in order to resum single-logarithmic terms and the type of double-sided singularities 
encountered in hadronic-jet profiles. 

The issues we have discussed are of some importance for QCD phenomenology. Before 
using any "safe" observable to test the theory or to measure as, one needs to identify the 
critical points of that observable and the expected behaviour in whatever order of pertur- 
bation theory is to be used. These points will need to be avoided in comparisons between 
fixed-order predictions and experiment. On the other hand, if resummed predictions can 
be obtained to single- logarithmic precision, then the behaviour at critical points can be 
used to provide interesting new tests of QCD, as was the case for resummation near the 
exclusive phase-space boundary ||||. 
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